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We introduce a new operator algebra for the description of the low-energy physics 
of one-dimensional, integrable, multicomponent quantum liquids. Considering the 
particular case of the Hubbard chain in a constant external magnetic field and with 
varying chemical potential, we show that at low energy its Bethe-ansatz solution can 
be interpreted in terms of the new pseudoparticle operator algebra. Our algebraic ap- 
proach provides a concise interpretation of and justification for several recent studies 
of low-energy excitations and transport which have been based on detailed analyses 
of specific Bethe-ansatz eigenfunctions and eigenenergies. A central point is that the 
exact ground state of the interacting many-electron problem is the non-interacting 
pseudoparticle ground state. Furthermore, in the pseudoparticle basis, the quan- 
tum problem becomes perturbative, i.e., the two-pseudoparticle forward-scattering 
vertices and amplitudes do not diverge, and one can define a many-pseudoparticle 
perturbation theory. We write the general quantum-liquid Hamiltonian in the new 
basis and show that the pseudoparticle-perturbation theory leads, in a natural way, 
to the generalized Landau-liquid approach. 
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I. INTRODUCTION 



For more than sixty years the "Bethe ansatz" (BA) [|T]|| has played a central role in the 
analytic solution of a variety of "integrable" Jl many-body problems in condensed matter 
physics and quantum field theory |l],|2|,|,|],^| . Initially applied to determine the ground 
state energies and spatial eigenfunctions, BA techniques were generalized to describe excited 
states, thermodynamics, and correlation functions and remain an active subject of 
study today. Spurred by a conjectured relationship to the microscopic mechanism for high- 
temperature superconductivity M , those particular integrable interacting electron models - 
such as the Hubbard Hamiltonian Q - which exhibit "Luttinger-liquid" behavior [0] have 
been extensively investigated [py9 ,1C| 11,|12,13 ,14, 15 . 



In particular, recent investigations have established the existence of Hubbard-model 
eigenstates possessing off-diagonal order (a feature essential for superconductivity) [|l(| and 
have shown that, although not themselves complete, the BA eigenstates can be extended 
to form a complete set of states for the one- dimensional Hubbard model ||17|| . Importantly, 
both these new results have been derived using powerful algebraic techniques, which rely on 
the symmetries of the Hubbard Hamiltonian. A vital and to date open issue is the extent to 
which these or similar algebraic techniques and symmetries can be used to describe specif- 
ically those excitations that dominate response and transport at low energies in integrable 
quantum liquids and related models. 



In the present article, which extends our earlier results ||18| , we show that there is a 
new operator algebra associated with the low-energy Hamiltonian eigenstates of integrable 
quantum liquids solvable by BA p],p|,^j8|JT3|] . We show that this algebra is expressed most 
naturally in terms of operators describing the "pseudoparticle" excitations introduced in 



several previous papers |Tj^^fIyJ^j2T|j2^| and proven to dominate low-energy transport 
and response functions. 

Whereas this previous work focused on the eigenstates and eigenfunctions of the pseu- 
doparticles, here we are able to work in terms of algebraic operators alone and hence to 



3 



obtain a much more compact and general representation. Considering the particular case 
of the Hubbard chain in a magnetic field and chemical potential, we demonstrate that the 
familiar BA solution can be interpreted naturally in terms of the pseudoparticle basis and 
that the Hamiltonian can be expressed simply in terms of operators in that basis. We then 
establish that a pseudoparticle perturbation theory, which is also naturally described in 
terms of the new operator algebra, can be used to study the low-energy excitations and the 
response and correlation functions. 

Our primary goal in the present article is to establish the utility of the pseudoparticle 
operator algebra for expressing the Hamiltonian and for calculating in perturbation theory 
the low-energy excitations of the one-dimensional Hubbard and related models. In a com- 



panion paper, which we shall henceforth call II ( p3|), we will apply the perturbation theory 
introduced here to the study of the Virasoro operator algebras of multicomponent integrable 
quantum liquids. Papers I and II contain a detailed exposition and extension of our earlier 
results 0. 

In Section II we introduce the Hubbard Hamiltonian in the presence of non-zero external 
magnetic field and chemical potential. We discuss prior results establishing that in this case 
the low-energy physics is dominated by a particular class of eigenstates: these lowest-weight 
states (LWS) of both the 77 spin and spin algebras ||,|l^,|4|,|5|,|6|,|^,|8|,|9|,|0| which refer 



to real rapidities [[l^]. We call these "LWS I" to distinguish them from the LWS associated 
with complex, non-real, rapidities, which we call "LWS II". In the sectors of parameter 
space with U(l) x U(l) symmetry the LWS II excitations have an energy gap and do not 
contribute to the low-energy physics. The LWS I can be described simply in terms of the 
pseudoparticles and the corresponding algebra. We demonstrate that in the pseudoparticle 
basis the ground state of the many-body problem is a "non- interacting" state, i.e., a simple 
Slater determinant of filled pseudoparticle levels. A central point is that for all canonical 
ensembles and symmetries with densities < n < 1 and spin densities < m < n the 
ground state is always an LWS I. We discuss the nature of the pseudoparticle Fermi sea 
and the low-energy excitations, showing that these excitations can be simply understood as 



pseudoparticle-pseudohole pairs and that they span the full low-energy Hilbert space. We 
introduce pseudomomentum-density operators and show that the operator expressions for 
the momentum, pseudoparticle number, and pseudoparticle-density fluctuations all have the 
expected forms. 

In Section III we relate the important rapidity numbers used in the BA to the 
pseudomomentum-distribution operators and use this relation to write the Hamiltonian 
(formally) in terms of the pseudoparticle operators. 

In Section IV we discuss the (potentially important) effects of normal-ordering and 
show that, in terms of the pseudoparticle operators, the normal-ordered Hamiltonian in 
the low-energy subspace has a simple form. This permits the introduction of a system- 
atic, non-singular pseudoparticle perturbation theory. Further, it establishes directly the 
universal form of the Hamiltonian for integrable multicomponent quantum liquids. These 
two results lead, in a natural way, to the one- dimensional Landau-liquid theory studied in 



Refs. KE , 151 rSLED ,EI , 22] . Our operator representation normal-ordered with respect to 



the ground state leads also to the results obtained from conformal- field theory [^,|T^], as 
we show in detail in paper II: the energy-momentum stress tensor and the generators of the 
Virasoro algebras can also be written simply in the pseudoparticle operator representation. 

In Section V we present a discussion and concluding remarks. We indicate how results 
derived for the Hubbard model can be extended to multicomponent BA solvable models 
|8| J13|JT8[| in which the theory has v > 2 independent branches of gapless elementary excita- 



tions. We contrast the perturbative character of the interaction effects in the pseudoparticle 
approach with the non-perturbative character of the usual electronic basis. We explain 
that the universal character of the present class of one- dimensional integrable quantum liq- 
uids can be understood in terms of a straightforward generalization to pseudoparticles of 



Wilson's renormalization group arguments ||31||: the pseudo- Fermi points of the pseudopar- 



ticles replace the particle Fermi surface and close to the pseudo-Fermi points only a few 
two-pseudoparticle scattering processes are relevant for the low-energy physics. Finally, we 
indicate briefly possible implications of our results for real quasi-one dimensional materials 



and for two-dimensional quantum systems. 



II. THE PSEUDOPARTICLE OPERATOR BASIS 



Consider the Hamiltonian describing the Hubbard chain in an external magnetic field H 



and with chemical potential \i [13,15,20,21 



H = -tJ2 [4 c i+i- + c )+i« c A + ^ EtsW " 1 / 2 }i4i c n - V2] - 2fifj z - 2^HS Z , (1) 



where t, U, fi, H, and /xq are the transfer integral, the onsite "Coulomb" interaction, the 
chemical potential, the magnetic field, and the Bohr magneton, respectively. The operator 
Cj a (cj a ) creates (annihilates) a spin a electron at the site j, and 
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(2) 



are the diagonal generators of the SU(2) rj spin and spin algebras, respectively p^6| , p!7| , |29| , |30 



As usual, a refers to the spin projections a =| , j, when used as an operator or function 
index and is given by o = ±1 otherwise. In Eqs. (2) 

N a = E c) a C ja = 4v C fe'<T = ^(0) , ^(k) = 4'+ka C k'a , (3) 

j k> k' 

is the number operator for spin a electrons. The operator p a (k) in Eq. (3) is the one-pair 
electron operator associated with the spin a fluctuations. Here c ka . (c^) is the momentum 
Fourier representation of the operator c] CT (cj CT ). 

The one-dimensional Hubbard Hamiltonian (1) describes an interacting quantum system 
of iVj_ down-spin electrons and iVf up-spin electrons on a chain of iV a sites with lattice 
constant a. Henceforth we employ units such that a = t = fio = h = l. Introducing 
n a = N a /N a and n = + n-f, we define {kp a — 7rn a ) and the total Fermi momentum by 
kp = im/2. 

In the absence of the chemical-potential and magnetic- field terms the Hamiltonian (1) 
has SO {A) = SU(2) ® SU(2)/Z 2 symmetry 0,0,^),gS],gB|,gS|,ggg . When N a is even, 
the operator f] z + S z (see Eq. (2)) has only integer eigenvalues and all half-odd integer 



representations of SU(2) <g>SU(2) are projected out jr^j28|] .The two SU(2) algebras - 77 spin 



and spin - have diagonal generators given by Eq. (2) and off-diagonal generators |TT|,|2"8 



v = Y,(-i) j mi, ^ = E(-i)MisV ( 4 ) 
3 3 

and 

S = 12 c h c n > ^ f = 12 c )i c 3i > ( 5 ) 
3 3 

respectively. 

In the presence of both the magnetic field and chemical potential terms, the symmetry 
is reduced to U(l) <8> U(l), with f) z and S z commuting with H. Hence at any fixed value of 
U, the parameter space of the Hamiltonian (1) is two dimensional [v dimensional for the 
multicomponent case when the number of subalgebras is v and the symmetry is [[/(l)]"). 
According to these values, the system has different symmetries as follows p9| : when rj z 7^ 
and S z ^ the symmetry is U(l) <g> U(l), for rj z = and S z it is SU(2) ® C/(l), when 
?7 2 7^ and ^ = it is Z7(l)(8)5'?7(2), and at r\ z = and 5^ = the symmetry is SO (4). Note 
that the ?7(1)®C/(1) symmetry sector always corresponds to two non-zero eigenvalues of the 
diagonal generators, whereas in sectors of higher symmetry, one or both of these eigenvalues 
vanish. 

Although the two conserved quantum numbers are usual taken (as in Eq.(3)) to be 
and iVj_, an alternative choice is N p and N az , i.e. N = N^ + N^ = /3 P (0) and N az = — N± = 
p CTz (0), where the charge and spin one-pair electron operators read 

p p {k)=p\k)+p l {k), p <T *{k)=p [ {k)-p l {k). (6) 

The eigenvalues r] z and S z can be expressed in terms of the conserved numbers, as 
shown by Eq.(2). This holds true also for the v eigenvalues of the diagonal generators of v 
multicomponent integrable systems. In these systems, in addition to the chemical potential 
and magnetic field, the system can include other external fields associated with the additional 
conserved quantum numbers. Note, however, that the number v of fields always equals the 
number of types of particles of the quantum liquid. 



In the above discussion we have worked in the parameter space specified by eigenvalues 
(Vz,S z ), corresponding to the canonical ensemble. Equivalently, we could choose the two 
external "fields" - H and \i - as parameters, or in fact use a "mixed" representation. We 
remark that in the papers p^|l5pO| , pTf a mixed representation involving the magnetic field 



H and the density n (i.e. number of particles N) was often used. 

Historically, most investigations of the Hubbard chain and related models have consid- 
ered the higher- symmetry sectors. For instance, in the case of the Hubbard model studies, 
have typically focused on the zero magnetic field (H = 0) case of symmetry £7(1) ® SU(2) 
(or 50 (4)). By considering the critical point of the one- dimensional Hubbard model in a 
magnetic field, Frahm and Korepin |13[ were able to discover interesting new features that 
remained hidden in previous studies of the model at zero magnetic field. For instance, in 
contrast to the H = case, where the spinon gapless excitations are LWS II, at H > (and 
n ^ 1) all gapless Hamiltonian eigenstates are LWS I and the form of the corresponding BA 
equations is much simpler. (The BA equations associated with the LWS II are in large num- 
ber and very complicated - see, for example, Refs. JT7||Z8| - whereas the BA equations which 



describe the LWS I are only two and have the same structure as the ground-state equations 
of Lieb and Wu m|I3|,|T3J.) However, as discussed in detail in Ref. ]nj, the zero magnetic-field 



limit of the H > low-energy expressions provides the correct H = exponents |13| and 



low-energy quantities [ 15| . Similarly, several papers focusing on the lowest symmetry region 



of parameter space have developed a pseudoparticle description and studied its relation to 



the Landau-liquid properties of the low-energy physics |[l^ , p!5] , p!8| , p^ , p0| , |^ , p^ , |23[ . 

We consider here the U(l) <S> U(l) sector which refers to electronic densities < n < 1 
and spin densities < m < n. This sector contains the "regular" BA states which are LWS 
of both algebras [[H]]. We concentrate our attention on the LWS I. In this U(l) ® U(l) 
sector 50(4) multiplets are generated by acting with the raising operators (4) and (5) on 



the corresponding LWS (for both LWS I and II). The total number of states is 4 a JT^ . 

It is important to reiterate that our focus on the U(l) <8> U(l) symmetry sector, in 
which the gapless BA Hamiltonian eigenstates are easier to describe, does not prevent us 
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from obtaining relevant information about the sectors of higher symmetry. This is because, 
as mentioned above, the limits of the expressions derived for the physical quantities in 
the lowest-symmetry parameter space generally approach the correct values for the same 
quantities in the sectors of higher symmetry ||15|| . When the transition from the sector of 
parameter space of lowest symmetry to the higher-symmetry sector does not involve the 
opening of a gap in one of the two elementary excitation branches, the values in the higher 
symmetry sector are obtained by taking in the expressions obtained in the U(l) <8> U(l) 
sector the limiting values for the eigenvalues of the diagonal generators which characterize 
the higher-symmetry sectors. For example, taking the limit S z — > in the expressions for 
the physical quantities derived in the present U(l) ® U(l) sector leads to the correct S z = 
expressions of the U(l) ® SU(2) (or 5*0(4)) sector. Although the metal-insulator transition 
can cause some subtleties pO| , |2T|j29[| in determining the n = 1 expressions - corresponding 
to the symmetry SU(2) <g> U(l) (or SO (A)) - these can also be obtained from the rj z — > 
values. 

As noted above, in the general v component case, in the sectors of lowest symmetry 
[^7(1)]^ the number v of U(l) subalgebras also provides the number of independent gapless 
branches of elementary excitations As we discuss below, each of these u branches is 



characterized by a different "color" quantum number, which we shall call a. These colors 
label the pseudoparticle operators and define "new" orthogonal directions in the low-energy 
Hilbert space. Furthermore, the numbers of pseudoparticles and pseudoholes are good quan- 
tum numbers. Therefore, at constant values of the particle-particle interaction U, we may 
use the v numbers of a pseudoparticles or the v numbers of a pseudoholes as alternatives to 
the v particle numbers N a associated with the operators (3), to define the parameter space. 

We return now to the specific case of the Hubbard model and study in detail the operator 
algebra introduced in Ref . Jl8| , which generates the Hilbert sub-space spanned by the regular 



BA eigenstates described by real rapidities (LWS I). 

A central point is that in canonical ensembles of the U(l) <S> U(l) sector both the non- 



LWS multiplets and the LWS II have energy gaps relative to the ground state fll5l , |18| , |29 



Therefore, for energies smaller than those gaps, the Hilbert subspace spanned by the regular 
BA LWS I coincides with the full Hilbert space of the quantum problem. The perturbation 
theory which we introduce in Sec. IV refers to that Hilbert space corresponding to energy 
scales smaller than the above gaps. We note, however, that LWS I contain both high- and 
low-energy states. In this section we study the entire Hilbert subspace spanned by the LWS 
I, but stress that it coincides with relevant, low-energy Hilbert space of the quantum problem 
only at energy scales smaller than the above gaps. 

The Hilbert subspace spanned by the LWS I can be generated by acting on the vacuum 
state with the pseudoparticle algebra we introduce below. The LWS I pseudoparticles are 
massless in the sense of field theory (describe gapless modes); of course, they have finite static 
and transport "condensed- matter" masses [[1^,|18[], which are given in II (for the transport 
masses see Eq. (71) of II.) For completeness we note that a more general Landau-liquid 



theory including the LWS II states can be constructed p2| . These LWS can also be described 
by pseudoparticles, but these new pseudoparticles have a mass in the field-theory sense (gap, 
in the present condensed-matter language). 

Although the pseudoparticles associated with the LWS I are the transport carriers at 
low energy [ 15] , they refer to purely non-dissipative excitations, i.e. the Hamiltonian com- 



mutes with the currents in the subspace spanned by the LWS I [|L5|,[^]. Therefore, these 
pseudoparticle currents give rise only to the coherent part of the conductivity spectrum, 
i.e. to the Drude peak |15l , |T8| , |23|| . The non-coherent part is associated with the massive 



pseudoparticles which describe the LWS II and some of the non-LWS multiplets |32| . In the 



Hilbert subspace spanned by those excitations, the Hamiltonian does not commute with the 
current operators. In addition, and to make life more difficult (and interesting), the c (s) 
excited LWS I states ( (18) below) disappear at n — 1 (or H — 0) and the corresponding 
massive pseudoparticle becomes massless (U(l) to SU(2) transition) p!5| , |32|j33| . Fortunately, 
for most of our discussion, this subtlety will be irrelevant. 

The new operator algebra studied in this paper generates all the LWS I from the electronic 
vacuum \ V). At constant values of U and S z < 0, the vacuum \ V) corresponds to the limit of 
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vanishing electronic density n — > 0. This is the S z < 0, U(l) (g> U(l) vacuum which refers to 
the limits — > and — ► with njn^ > 1. We note that when this vacuum is obtained 
from the ground state (given in Eq. (11) below) by taking that limit, in the case of some 
matrix elements this does not commute with the limit n-f — > and — > with njn^ = 1, 
which defines the = 0, £7(1) <8> SU(2), vacuum. In turn, these two limits lead to the same 
S z = vacuum when they refer to the ground states alone (see Eq. (11) below), i.e. not to 
the limiting value of matrix elements involving these states. This is revealed by the study 
of the electron - pseudoparticle canonical transformation which will be presented elsewhere 







The new operator algebra involves two types of pseudoparticle creation (annihilation) 
operators b qa (b qa ) which obey the usual fermionic algebra fL8| 



Here a refers to the two (or v, in the general case) pseudoparticle colors c and s |15,|l8,29 . 
The discrete pseudomomentum values are 

where J" are consecutive integers or half integers. (In the case of the Hubbard chain we 
have that = Ij and J| = J@ [with (3 — j], where Ij and Jp refers to the notation of Lieb 
and Wu [f|.) There are N* values of Jj*, i.e. j — 1, iV*. An LWS I is specified by the 
distribution of N a occupied values, which we call a pseudoparticles, over the iV* available 
values. There are N* — N a corresponding empty values, which we call a pseudoholes. We 
emphasize that the BA wave functions vanish for double-occupied configurations of the 7° 
quantum numbers ||. 

In the present U(l) ® U(l) sector of parameter space, the numbers of a pseudoholes 
iV* — N a , such that A*"* > N a , are determined by the corresponding eigenvalues of the two 
(or v) diagonal generators. In the case of the Hubbard chain we have that 

-V, = \[K - N c ] , - Sz = ^[N:-N s ], (9) 
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where 



N* c =N a ; N C = N; N* s = iV T ; N s = N± . (10) 

The numbers ij are integers (or half integers) for N s even (or odd), and IJ are integers 
(or half integers) for N* odd (or even) ||. Since only single and zero occupancy of the 
values IJ are allowed, only pseudoparticles of the color a can occupy the states labeled by 
the numbers IJ. Therefore, the pseudoparticles have a fermionic character, as assured by 
the anticommuting algebra (7). For the Hubbard model the BA spatial wave function for the 
LWS I depends on the quantum numbers IJ through two sets of real numbers, which many 
authors call rapidities. Note that for the LWS II these rapidities are complex, non-real, 
numbers. The expression of the spatial wave functions in terms of the quantum numbers IJ 
requires the solution of two systems of algebraic equations which define the two rapidities 
as functions of the quantum numbers I J PJ^ Jl3|j2]| . Although the expression of the spatial 



wave function for the LWS I in terms of the quantum numbers IJ requires the solution of 
the above systems of equations, which constitutes a problem of considerable complexity, 
the description of these eigenstates in the basis of the above BA quantum numbers IJ that 
diagonalize the quantum liquid is much simpler. Historically, the Hamiltonian eigenstates 
were introduced in terms of the spatial BA wave functions . The diagonalization of the 
problem leads then to the two systems of algebraic equations. In the case of LWS I these 
equations introduce the integer or half-integer quantum numbers IJ which describe these 
eigenstates. 

One of the principal advantages of the algebraic approach we use in this paper is that 
it permits the description of the LWS I in terms of the quantum numbers IJ and does 
not require the spatial wave-function representation. In the basis associated with these 
quantum numbers the description of the LWS I does not involve the rapidity numbers. As 
in the case of the spatial wave functions, the expression of the energy in terms of the quantum 
numbers IJ involves the rapidities. In Sec. Ill and Appendix A we will consider the above 
algebraic equations for the case of the LWS I of the Hubbard chain. The introduction of 
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a suitable operator representation reveals that the rapidities are the eigenvalues of rapidity 
operators which in the pseudoparticle basis determine the Hamiltonian interaction many- 
pseudoparticle terms. 

The above electronic vacuum state \V) has no electrons. At fixed values of the onsite 
interaction U, the ground state of a canonical ensemble characterized by eigenvalues r\ z < 



and S z < can be constructed from that vacuum. Following Refs. |I^ , p!8|j20| , |21|j29|| and 



as shown elsewhere |30|j32|], it corresponds to filling symmetrically around the origin N a 
consecutive 7™ values of all colors a. In addition, the ground state of canonical ensembles of 
all symmetries and such that < n < 1 and < m < n are always LWS I. This holds even 
at the SO (4) "point", where n = 1 and m = 0, and the available Hilbert space contains only 
one LWS I, which is the corresponding ground state; in this case, all LWS excited states 
are LWS II, as we discuss below. The ground state associated with a canonical ensemble of 
(r] z , S z ) values has the form 

io;th,s,)= n [ n *u\v), (ii) 

where when N a is odd (even) and 7? are integers (half integers) the pseudo-Fermi points 
are symmetric and given by 

= -«£ = £[*a " 1] ■ (12) 

On the other hand, when N a is odd (even) and 7j* are half integers (integers) we have that 

q% = = *]> (is) 

or 

$2 = -£-[Na-2], -qti = ^-N a . (14) 

Similar expressions are obtained for the pseudo-Brioullin zones limits if we replace in 
Eqs. (12) - (14) N a by N*. 
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The simple form of the ground-state expression (11) has a deep physical meaning. It 
reveals that in the pseudoparticle basis the ground state of the many-electron quantum 
problem is a "non-interacting" pseudoparticle ground state of simple Slater-determinant 
form. However, that the numbers I" of the RHS of Eq. (8) can be integers or half integers for 
different ground states which, for example, differ by 1 in one of the eigenvalues of the diagonal 
generators (9), makes the problem much more involved than a simple non- interacting case. 
This change in the integer or half-integer character of some of the numbers J" of two ground 
states, shifts all the occupied pseudomomenta and leads to the orthogonal catastrophes 
0. This is the reason for the absence of quasiparticle peaks in the single-particle spectral 
function. On the other hand, if we consider the evaluation of quantities which do not change 
the eigenvalues of the diagonal generators (9), the integer or half-integer character of the 
numbers Jj* is, in the thermodynamic limit, irrelevant. 

In some problems it is important to distinguish the two situations = —qpa and 
#Fa 7^ ~Qfcx- For instance, the total- momentum expression (22) below reveals that when 
Qfo 7^ — QfI ^ ne a pseudoparticles give a finite contribution to the momentum of the ground 
state (11). However, in many evaluations we can assume the simplest case qp^ = —qp2 and 
similarly for the pseudo-Brillouin zones limits. This is because in either case qp^ — qp~J = 
2qFa — 27r /N a and q^ — q^ = 2q a — 2n/N a , and, in addition, except for terms of order 
1/N a , we have that qp^ = -qp2 = QFa and q^ = -gH = g Q) where 

' • (15) 



The representation (8) was used by Yang and Yang in Ref. ||35|| . In the thermodynamic 



limit we take the continuum limit qj — > q and the set of possible pseudomomenta (8) maps 
into a continuum domain of q values. For the calculation of many quantities we can then re- 
place ±<Zi^ and ±gi ± ' ) by the values qp a and q a , respectively, given in Eq. (15). For instance, 



in Refs. |]T5|J20| , |2T|| that approximation was used in the evaluation of the Landau-functional 
expansions because this leads to the correct values for the corresponding pseudoparticle 
bands and / functions. In contrast, however, when we insert in that functional particular 
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forms of the pseudoparticle deviations which refer to changes in the numbers N a , the integer 



or half-integer character of the quantum numbers 7? plays a crucial role . 

As in the case of Landau's Fermi liquid theory p6| , |37| , we will use the exact ground state 
(11) of the quantum problem as reference state. In the pseudoparticle basis we will write 
all operators in normal order relative to that ground state. The corresponding universal 
form of the normal-ordered quantum-liquid Hamiltonian will confirm that in that basis the 
present quantum system is a Landau liquid []l4j , |l^[l9| , ^|21 1 . 

It is useful to know how many LWS I and LWS II there are for given numbers r] z and S z . 
The total number of LWS, i.e. of regular BA Hamiltonian eigenstates, was evaluated in Ref. 



17 1 . The following numbers refer both to canonical ensembles of the present 7/(1) ® U(l) 
sector and of the U(l) <g> SU(2), SU{2) ® 17(1), and SO (A) sectors when < n < 1 and 
< m < n. Only in the 77(1) (g> U(l) sector do all the LWS II states have an energy gap 
relative to the ground state (11). The number of LWS I and LWS II is given by 



and 



(16) 



N a 




n; 


n s + n: 




N s 



N a 
N* 



N a 



N a 
TV. -2 



N a 
N* + 1 



N* 



N s 



N n 



N s + N* 



n; 

N s 
(17) 



respectively [^]. The square brackets in the above equation refer to the usual combinatoric 
coefficients. Adding the numbers (16) and (17), one recovers the total number of regular 
BA Hamiltonian eigenstates for given values r] z and S z JT7J- 

Since the colors a and the pseudomomentum q are the only quantum numbers involved 
in the description of the pseudoparticles whose distributions define the LWS I, these dis- 
tributions can be generated by applying to the vacuum \V) the J2 a N a creation operators 
b^a- The resulting Hamiltonian eigenstates, which are LWS I and of total number given by 
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(16), include the ground state (11) and the excited states which can be generated from it 
by pseudoparticle-pseudoholes processes and have the form 



pfl 



1»S*)= II til 6L6«J|0;%,^>, (18) 



1i 

a=c,s i,j=l 



where qj (g^) defines the different locations of the pseudoparticles (pseudoholes) relative to 
the reference state (11) and N^ h is the number of a pseudoparticle-pseudohole processes. 
5*0(4) multiplets are generated by acting the raising generators (4) and (5) of the rj spin 
and spin algebras, respectively, onto the LWS I (11) and (18). As is shown in Ref. |30|, in 



contrast to the higher-symmetry sectors, in the present U(l) ® C7 (1) sector we have that 
non-LWS multiplets with values of r] z and S z as in the corresponding canonical ensemble 
have an energy gap relative to the ground state (11). 

The LWS I states of the form (11) and (18) and of total number given by (16) constitute 
a complete orthonormal basis which spans an important Hilbert subspace, which we call Tij. 
At energy scales smaller than the gaps for non-LWS multiplets |JD| and LWS II, TCj repre- 



sents the full accessible Hilbert space. We emphasize that the generators which transform 
the ground state (11) into the excited eigenstates (18) are products of one-pseudoparticle 
operators of the form b q+ka b qa , where < \q\ < qp^ and qp^ < \q + k\ < 

Equations (11) and (18) reveal that in the basis associated with the pseudoparticle op- 
erators b' and b qa , all Hamiltonian eigenstates of the many-electron system which are LWS 
I have a "non-interacting" form. Since for a given canonical ensemble of values (r] z ,S z ) all 
these excitations can be generated by successive applications of the operator b q+ka b qa to 
the ground state (11), the expression for this operator in the usual electronic basis would 
provide interesting information about the nature of the states (18) from the point of view 
of electronic configurations: while at large ±2qp a momenta and low energy the one-pair a 
pseudoparticle-pseudohole Hamiltonian eigenstates can correspond to multipair electronic 
excitations |22|,|34|], at small values of the momentum the operators b q+ka b qa are a superpo- 
sition of the two one-pair electronic (or particle) operators (4) P9|,51|. Fortunately, the BA 



solutions are naturally expressed in terms of the pseudoparticle basis. Therefore, although 
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we do not know, in general, the precise form of the above generators in the electronic rep- 
resentation, these solutions provide the expression for the Hamiltonian and other operators 
in the pseudoparticle basis, and we can extract relevant information about the quantum 
system without describing the problem in terms of electronic configurations. An impor- 
tant limitation, however, is that the BA solution provides the pseudoparticle expressions for 
some operators only. The general problem of the electron - pseudoparticle canonical trans- 
formation will be studied elsewhere p9U34). The fact that the BA solutions are not most 
naturally expressed in terms of the original electronic basis is the main reason why it has 
been difficult, in previous work, to extract the information concerning correlation functions 
and matrix elements contained in the BA solutions. 

Despite the non-interacting form of the Hamiltonian eigenstates (11) and (18), we find in 
Sec. IV that the normal-ordered Hamiltonian includes pseudoparticle interaction terms and 
is, therefore, a many-pseudoparticle operator. However, these pseudoparticle interactions 
have a pure forward-scattering, zero-momentum transfer, character. This agrees with the 
Landau- liquid studies of Refs. |r^ , |T5|]20|pi| , which established this result using eigenenergies 



and eigenfunctions rather than operators. 

The expression of the Hamiltonian (1) in Tij involves exclusively the two (a = c,s) 
pseudomomentum distribution operators 

N a (q) = bl a b qa , (19) 

which play a key role in the present basis. The operators (19) commute with each other, 
i.e. [N a (q), N a /(q')} = 0. Note, however, that the two (or v) pseudomomentum distribu- 
tion operators (19) are not the two momentum distribution electron (or particle) operators 
N a (k) = c\ a Cka- The operator N a (q) (N a (k)), which has a simple form in the pseudoparticle 
(electronic) basis, has an involved expression in the electronic (pseudoparticle) basis p3fl . 

Since in the pseudoparticle basis the Hamiltonian expression involves only the operators 
(19), it follows that in Tij the Hamiltonian commutes with these operators. These play a 
central role in this Hilbert subspace because all the Hamiltonian eigenstates which are LWS 
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I are also eigenstates of N a (q). These LWS I are of form (11) or (18) and, therefore, obey 
eigenvalue equations of the form 

N a (q)\Vz,S z ) = N a (q)\ Vz ,S z ), (20) 

where N a (q) represents the real eigenvalues of the operators (19), which are given by 1 and 
for occupied and empty values of q, respectively 

It follows from Eq. (9) that in 7ix and in the pseudoparticle basis the diagonal generators 
(2) of the two U(l) diagonal subalgebras become 

V* = -~ EI 1 - N c (q)\ , S z = -\ Y\l ~ N s (q)} . (21) 

Z q Z q 

When we take the zero limit in one or several of the eigenvalues of the diagonal generators, 
i.e. limit 7] z — > or (and) S z — > 0, we approach a boundary line sector of parameter space 
of higher symmetry. Following Eqs. (9) and (21), this corresponds to taking one or several 
of the limits N a — > N*. If, for example, we take this limit for one of the a pseudoparticle 
branches only, since the corresponding number of a pseudoholes vanishes, Eqs. (15) and 
(16) reveal that the number of a configurations reduces to 1. This is a filled sea of a 
pseudoparticles similar to that of the corresponding ground state (see Eqs. (11) and (15)). 
In particular, when we take the limit N a — > N* for all a pseudoparticle branches, Eq.(16) 
gives one single configuration which is the 5*0(4) half-filling and zero-magnetic field ground 
state [p9fl . In this case there are no excited LWS I (18) and, except for the ground state, 
all LWS are LWS II p8| , p9| . On the other hand, it also follows from Eq. (16) that taking 
the limit N a — >• in all pseudoparticle branches corresponds also to a single configuration. 
In this case all a bands are empty, i.e. full of pseudoholes. This single state is the vacuum 
\V), which is both the electronic and pseudoparticle vacuum. 

These arguments show that, when we reach a parameter-space sector of higher symmetry 
and/or a new phase by changing the eigenvalue of one of the diagonal generators (21), the 
corresponding a branch of pseudoparticle-pseudohole eigenstates disappears. In the case 
of a phase transition, the collapsing branch is not replaced by a new gapless branch of 
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LWS II, whereas in the "microscopic transitions" |15|j29|] one of the branches of LWS II 
becomes gapless. In the latter case the number of gapless branches remains the same, but 
the symmetry is higher. In the sector of parameter space of lowest symmetry we have that 
the number of a configurations in (16) is larger than one for all colors a. Moreover, all 
non-LWS multiplets |JIJ an d LWS II, the study of which we have omitted from this paper, 
have an energy gap []l5 , |20 ,29l,32 



Equation (21) gives the expression of the diagonal generators (2) in the pseudoparticle 
basis. Combining the operator representation that we have just introduced with the infor- 
mation contained in the BA solution, it is straightforward to write other simple operators 
in the pseudoparticle basis. The momentum P and and the number of a pseudoparticle 
operator N a , for example, read 

P = £giV a (g), N a = Y,N a (q), (22) 

q,a q 

respectively. These expressions were derived by combining Eqs. (10) and (20) with Eqs. 



(8) — (10) of Ref. f2"0fl . Note that these operators contain non-interacting pseudoparticle 
terms only. 

To close this section, and in order to clarify the relation of our operator representation 



to the Landau-liquid functional of Refs. [Tj.[T^.|2yj2T|] , we discuss the eigenvalues N a (q) of 
Eq. (20). From Eq. (22) we observe that these eigenvalues obey the following normalization 
equation 



N a = = / dqN a (q) . (23) 

q 27T Jq { a > 

For instance, in the case of the ground state (11) we have that 

N° a (q) = 6(#J ~q), 0<q<q^ 

= ®(q-q ( il), q { - ) <q<0, (24) 

where q^2 is defined in Eqs. (12) — (14) and q^ is given by similar equations with N a 
replaced by N*. 
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Using the ground state (11) as the reference state, among all the excited eigenstates 
of the form (18) we will be mostly interested in low-energy excitations. These involve a 
redistribution of a small density of pseudoparticles relative to the distribution (24). Let us 
introduce the normal-ordered pseudomomentum distribution operator 

: N a (q) := N a (q) - N° a (q) , (25) 

which is such that (0;r] z ,S z \ : N a (q) : \0;r] z ,S z ) = 0. The normal-ordered operators (25) 
obey the eigenvalue equations 

:N a (q): \rj z , S z ) = 5N a (q)\ Vz , S z ) , (26) 

where \rj z ,S z ) denotes any Hamiltonian eigenstate of form (11) or (18). Here SN a (q) = 
N a (q) — N®(q) and, therefore, 5N a (q) = when \r] z , S z ) = |0; r/ z , S z ). The eigenvalues N a (q) 
and 5N a (q) are nothing but the pseudomomentum distributions and deviations, respectively, 
of the Landau-liquid theory studied in Refs. JT^| , |T5] , [r9] , pD] , pT|] . Equations (25) and (26) imply 
that these are the expectation values: 

N a (q) = (r ]z ,S z \N a (q)\ Vz ,S z ), (27) 

and 

SN a (q) = ( Vz ,S z \:N a (q):\r ]z ,S z ), (28) 

respectively. Therefore, the introduction of the operator algebra (7), (11), and (18) clarifies 
the deep reasons for the validity of the one-dimensional Landau-liquid theory, as discussed 
in more detail in Sec. IV. 



III. RAPIDITY OPERATORS 

In this section we continue to restrict our study to Tij- As discussed above, in the case of 
LWS I, the BA solutions of the Hubbard model lead to two systems of algebraic equations. 
Given the configuration of quantum numbers that describes each Hamiltonian eigenstate of 
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form (11) or (18) (LWS I), these equations fully define the two sets of rapidities which deter- 
mine the corresponding spatial wave function. As in the case of the quantum numbers J", 
each of these types of rapidities is associated with one of the colors a. In the thermodynamic 
limit, we can take the pseudomomentum continuum limit qj —>■ q and the rapidities become 
functions of q which we call R a (q). If we combine the operator representation introduced 
in Sec. II with the properties of the BA solution, it is straightforward to show from the 
relation between the rapidities R a (q) and the pseudomomentum distributions (19) and (25) 
that these rapidity functions are nothing but real eigenvalues of rapidity operators R a {q) 
such that 

Ra(q)\Vz,S z ) = R a (q)\Vz,S z ). (29) 

The two rapidity operators R a (q) contain all information about the many-pseudoparticle 
interactions of the quantum-liquid Hamiltonian. There are two fundamental properties 
which imply the central role that the rapidity operators of Eq. (29) have in the present 
quantum problem: 

(a) Each of the normal-ordered rapidity operator : R a (q) : can be written exclusively in 
terms of the two pseudomomentum distribution operators (25); 

(b) The normal-ordered Hamiltonian can be written, exclusively, in terms of the two pseu- 
domomentum distribution operators (25), but all the corresponding many-pseudoparticle 
interaction terms can be written in terms of the rapidity operators : R a (q) :. It follows that 
in the Tij the rapidity operators commute with the Hamiltonian. 

In the thermodynamic limit the BA algebraic equations for the LWS I are replaced by a 
system of two coupled integral equations The standard treatment of the BA solutions 
is very lengthy, e.g., for each eigenstate we have to rewrite a new set of equations. This is 
because the eigenvalues R a (q) of the rapidity operators R a (q) are different for each state 
\Vz,S z ) (18). 

However, within the present operator description, we can introduce a single set of two 
general operator equations which apply to all eigenstates. Each of these equations defines 
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one of the rapidity operators R a (q) in terms of a q and a summation containing functionals of 
the rapidity operators and pseudomomentum-distribution operators. This system of coupled 
equations has a unique solution which defines the expressions for the rapidity operators 
in terms of the pseudomomentum-distribution operators. Although the structure of these 
equations leads to many universal features which we will discuss later, they involve the 
spectral parameters which are not universal but are specific to each model. Also the formal 
expression for the Hamiltonian in terms of the rapidity operators is not universal. However, 
the corresponding normal-ordered expressions do have a universal form for the integrable 
multicomponent systems, as we shall discuss in Sec. IV. 

For consistency with the previous notation Hl5|,^0|j2~lj , we use for the Hubbard chain 
the notation R c {q) = K{q) and R s {q) = S(q). The spectral parameters are the numbers 
(4:t/U) sin[K(q)] and S(q), which appear in Appendix A and in Eqs. (30) — (32) below in 
operator form. 

We have mentioned several times that the BA solution is most naturally expressed in the 
pseudoparticle basis. One reflection of this is the simple expression for the Hamiltonian (1) 



in that basis. Combining Eq. (29) with the energy expression (4) of Ref. which refers 
to the LWS I, the result is 

H = £ JV c (g){-2*cos[£(g)] - U/2} - 2^ z - 2^HS Z , (30) 

where the expressions of the diagonal generators are given in Eq. (21). This is the exact 
expression of the Hamiltonian (1) in Tij. At energy scales smaller than the gaps for the 
non-LWS multiplets and LWS II, (30) gives the exact expression of that Hamiltonian in the 
full Hilbert space. Despite its simple appearance, the Hamiltonian (30) describes a many- 
pseudoparticle problem. The reason is that the expression of the rapidity operator K(q) in 
terms of the operators N a (q) contains many-pseudoparticle interacting terms. 

To specify the Hamiltonian (30) completely, we must indicate the operator equations 
that define the rapidity operators in terms of the pseudoparticle momentum distribution 
operators. In the case of the Hubbard chain these two operator equations read 
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{K(q) - ^-E^(^)tan- 1 (S(g') - (At/U) sm[K(q)})]\ Vz , S z ) = q\ Vz ,S z ) (31) 



N a ql 



and 



A E iV c ( g ')tan- 1 (5(g) - (At/U) sin[i%')]) 



N a ql 



£ ^W) tan" 1 (i (5(g) - S(q')))}\ Vz , S z ) = q\ Vz , S z ) . (32) 



3' 



These equations fully define the rapidity operators in terms of the pseudomomentum distri- 
bution operators (19) and (25). 

Simple generalizations of these equations hold for other multicomponent quantum liquids: 
the many-pseudoparticle terms of the Hamiltonian can be expressed, exclusively, in terms 
of the rapidity operators. Moreover, the BA solutions always provide: (a) the expressions 
of the Hamiltonian in terms of the rapidity operators; and (b) v operator equations which 
define (implicitly) the expressions of the rapidity operators in terms of the pseudomomentum 
distribution operators (19). 

The forms of the eigenfunction equations (20) and (29) show that the operator equations 
(31) and (32) are equivalent to the BA equations presented previously |T5|J2y,|2T]| . Formally, 



these are obtained from Eqs. (31) and (32) by replacing the operators by the correspond- 
ing eigenvalues and projecting onto the various \r] z , S z ) states. As in the case of Eqs. (31) 
and (32), these equations are general and apply to all eigenstates. Their solution gives 
the rapidity-Landau functionals in terms of the pseudomomentum distributions (19). In- 
sertion into these functionals of the pseudomomentum distributions of a given Hamiltonian 
eigenvalue \r] z , S z ) leads to the corresponding rapidity eigenvalue of the RHS of Eq. (29). 

Unsurprisingly, it is difficult to solve the BA operator equations (31) — (32) directly and 
to obtain the explicit expression for the rapidity operators in terms of the pseudomomentum 
distribution operators (19). In contrast, it is easier to calculate their normal-ordered expres- 
sion in terms of the normal-ordered operators (25). In the ensuing section, we introduce the 
pseudoparticle perturbation theory which leads to the normal-ordered expressions for the 
Hamiltonian and rapidity operators. 
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IV. PSEUDOPARTICLE PERTURBATION THEORY 



In the pseudoparticle basis the normal-ordered rapidity operators : R a (q) : contain an 
infinite number of terms, as we shall demonstrate below. The first of these terms is linear 
in the pseudomomentum distribution operator : N a (q) : (25), whereas the remaining terms 

consist of products of two, three, , until infinity, of these operators. The number of 

: N a (q) : operators which appears in these products equals the order of the scattering in the 
corresponding rapidity term. 

A remarkable property is that in the pseudoparticle basis the seemingly "non- 
perturbative" quantum liquids become perturbative: while excitations associated with adding 
or removing of one electron have zero life-time and decay into collective pseudoparticle ex- 
citations, all low-energy Hamiltonian eigenstates of the system with one electron more or 
less are pseudoparticle-pseudohole states of the form (11) or (18). Furthermore, while the 
two-electron forward scattering amplitudes and vertices diverge, the two-pseudoparticle / 
functions (given by Eq. (42) below) and the corresponding two-pseudoparticle forward- 
scattering amplitudes, which were calculated in Ref. ||15|| , are finite. By "perturbative" we 



also mean here the following: since at each point of parameter space (canonical ensem- 
ble) the excited low-energy eigenstates are of form (18) and correspond to quantum- number 
configurations involving a density of excited pseudoparticles relative to the ground-state con- 
figuration (11), (24), the expectation values of the Hamiltonian in these states are functions 
of the density of excited pseudoparticles. This density is given by 

n ex = J2 n ex > ( 33 ) 

a 

where 

<*=^£[l-^a(?)]^(<?) (34) 

defines the density of excited a pseudoparticles (and a pseudoholes) associated with the 
Hamiltonian eigenstate \r) z ,S z ). 
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When all the densities n^ x are small, we can expand the expectation values in these den- 
sities. The perturbative character of the quantum liquid rests on the fact that the evaluation 
of the expectation values up to the n th order in the densities (33) requires considering only 
the corresponding operator terms of scattering orders less than or equal to n. This follows 
from the linearity of the density of excited a pseudoparticles, which are the elementary 
"particles" of the quantum liquid, in SN a (q) = (r] z ,S z \ : N a (q) : \r) z ,S z ) and from the form 
of (34). The perturbative character of the quantum liquid implies, for example, that, to 
second order in the density of excited pseudoparticles, the energy involves only one- and 
two-pseudoparticle Hamiltonian terms [|18|,^| , as in the case of the quasiparticle terms of a 
Fermi-liquid energy functional [p6|,p7| . 

We note parenthetically that Eqs. (33) and (34) do not apply to the excitations (A) 



studied in [[HJ and in II. In the case of these excitations, which are associated with changes 
in the numbers N a , the densities (33) and (34) are replaced by the density of "removed" or 
"added" pseudoparticles. In this case the change in the integer or half-integer character of 
the quantum numbers must be taken into account. 

The eigenvalue equations (26) imply that the problem of using the rapidity Eqs. (31) 
and (32) to derive the expression of the operators : R a (q) : in terms of the operators 
: N a (q) : is equivalent to the problem of evaluating the corresponding expansion of the 
rapidity eigenvalues 5R a (q) in terms of the pseudoparticle deviations SN a (q) (28). This last 
problem, which leads to the Landau-liquid expansions, was studied in Ref. |^TJ for the case of 
the Hubbard chain. Furthermore, we emphasize that it is the perturbative character of the 
pseudoparticle basis which justifies the validity of these Landau-liquid deviation expansions 
0,110,11. 

Based on the connection between the two problems we can derive the expressions for the 
rapidity operators : R a (q) :. This corresponds to expanding the expressions of the operators 
: R a (q) '■ in terms of increasing pseudoparticle scattering order. It is convenient to define 
these expressions through the operators : Q a (q) :. These are related to the rapidity operators 
as follows 
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:R a (q):=R a (q+:Q a (q):)-R a (q), (35) 

where R Q a (q) is the ground-state eigenvalue of R a (q), i.e. 

R a (q)\0; S z ) = R° a (q)\0; Vz , S z ) . (36) 

The operators : Q a {q) '■ contain the same information as the rapidity operators, and involve, 
exclusively, the two-pseudoparticle phase shifts $ aa '(q, q') defined in Ref. |^T[. Following the 
related studies of Refs. |T5| , pO| , |2T|] , in Appendix A we introduce Eq. (35) in the BA equations 
(31) and (32) and expand in the scattering order. This leads to 

oo 

:&.(*) :=£<K } (ff), (37) 

i=l 

where i gives the scattering order of the operator term Q$(q). For example, for the first- 
order term we find 

Q { « ) (q) = ^Y,^ a '(q,q'):N a ,(q'): . (38) 

We emphasize that while the expressions for the phase shifts $ aa ' (q, q') are specific to each 
model because they involve the spectral parameters |[21|| , the form of the operator term 



Qa\q) (38) is universal and refers to the multicomponent quantum liquids discussed in 
Sees. I and V. 

All the remaining higher-order operator terms of expression (37), Q^(q), can be obtained 
from the rapidity equations (31) and (32). For simplicity, we provide only the expression 
of the first-order operators (38) and of the second-order operators, which are given in Eqs. 
(A20)-(A23) of Appendix A. In that Appendix we evaluate the expressions of these operators 
and use the Hamiltonian expression (30) in terms of the rapidity operators to derive the 
following expression for the normal-ordered Hamiltonian : 

oo 

;H;=Y,H {i) , (39) 
i=l 

where, to second pseudoparticle scattering order 
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H (1) = £e Q (g):A> a (g): 



q.a 



H {2) = ^EE ^aa'(g,?') : iV Q (g) :: A> Q ,(g') : . (40) 



Here (40) are the Hamiltonian terms which are relevant at low energy |18|,^3j. Furthermore, 



it is shown in Refs. |18| , |23| that at low energy and small momentum the only relevant term 
is the non-interacting term This property justifies to the Landau-liquid character of 

these systems and plays a key role in the symmetries of the critical point |1S|]23|| . 

The form of the normal-ordered Hamiltonian (39) — (40) is universal for the integrable 
multicomponent quantum liquids. On the other hand, the expressions for the pseudoparticle 
bands e a (q) involve the spectral parameters and are specific to each model. For the case 
of the Hubbard chain they are defined in Ref. p0| . The c and s pseudoparticle bands are 
shown in Figs. 7 and 8, respectively, of that reference. 

All Q^(q) terms of the RHS of Eq. (37) are such that both the / functions of the RHS 
of Eq. (40) and all the remaining higher order coefficients associated with the operators 
jjW f order i > 1 have universal forms in terms of the two-pseudoparticle phase shifts and 
pseudomomentum derivatives of the bands and coefficients of order < i. This follows from 
the fact that the S'-matrix for i-pseudoparticle scattering factorizes into two-pseudoparticle 



scattering matrices, as in the case of the usual BA S'-matrix p^j38|j39|j40[1 . For example, 



although the second-order term H^> of Eq. (40) involves an integral over the second-order 
function Q^\q) (see Eq. (A24) of Appendix A), this function is such that can be 
written exclusively in terms of the first-order functions (38), 

H {2) = ^MQSKq) : N a (g) : + ^£^ £ [ffOWP. (41) 

q,a Z7r a * j=±l 

as shown in Appendix A, and the four (or, in the general case, v x v) "Landau" / functions, 
faa' (?) <?') 5 have universal forms which read 

faa'(q, q') = 2nv a (q)® aa <(q, q') + 2-KV a ,{q') $ a >a{q', q) 

+ £ £ 2Trv a/l ® a/la (jq Fa , l ,q)$ a „ a ,(jq Fa „,q f ) , (42) 

3=il a"=c,s 
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where the pseudoparticle group velocities are given by 

«-(?) = -^-. ( 43 ) 
and depend on the spectral parameters. In particular, the velocities 

v a = v a (q Fa ) , (44) 

play a determining role at the critical point, representing the "light" velocities which appear 
in the conformal-invariant expressions |ji~3l, 18|,P3| . In the case of the Hubbard chain, the 



velocities (44) are plotted in Fig. 9 of Ref . p0| . 

We note that the Hamiltonian term in (40) has, from the point of view of the 
pseudoparticle basis, a non-interacting character. However, in the electron basis, is of 
interacting character, as revealed by the U dependence of the bands e a (q) p0| . Furthermore, 
as we mentioned above, both the two-pseudoparticle / functions (42) and forward-scattering 
amplitudes ]I5] are finite. This is in contrast to the non-perturbative electronic basis, where 
the two-electron forward-scattering vertices and amplitudes diverge. 

Obviously, Eqs. (27) and (28) imply that the operator expressions (35) and (37) — (40) are 
fully equivalent to the corresponding Landau expansions already studied in Refs. J11^ . |2T)| . |2"1"[ 
for the case of the Hubbard chain. Thus we have achieved one of the aims of the present 
paper, which is to use the pseudoparticle operator algebra to justify the validity of the 
Landau-liquid properties of the quantum problem already studied in these earlier papers. 

We emphasize that in the present U(l) <g> U(l) sector of parameter space and at energy 
scales smaller than the gaps for the non-LWS and LWS II, Eqs. (30) and (39) — (40) refer 
to the expression of the full quantum-liquid Hamiltonian. In the electronic basis this is 
given by Eq. (1). Our operator representation leads, in a natural way, to the low-energy 
spectrum studied in Refs. [|T|,[]^,|1]|.^.|T]] . In particular, the study of Refs. p^JT5|j20|j2T 



shows that the low-energy of the Hamiltonian (39) — (40) is, in many points, similar to that 
of a Fermi liquid, with the pseudoparticles playing the role of the quasiparticles and the 
colors a of the spin projections a. For example, the static susceptibilities can be calculated 
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as in a Fermi liquid Plf , and involve, exclusively, the velocities (44) and the expressions of 



the two-pseudoparticle phase shifts at the pseudo-Fermi points. Also the low-temperature 
thermodynamics can be studied as in a Fermi liquid, being described by pseudoparticle 
Fermi-Dirac distributions of the form [EU[ 



N ^= 1 + e X)/ kB T - («) 

As in a Fermi liquid, the low-temperature specific heat is linear in the temperature and 
involves the static masses m* = where the pseudo- Fermi pseudomomenta and veloci- 
ties are defined in Eqs. (15) and (44), respectively. The low-frequency and low-momentum 
dynamical properties can be studied by means of kinetic equations [0,0, again as in a 
Fermi liquid. This allows the evaluation of the pseudoparticles elementary currents and 
transport masses. The latter masses, the expressions for which differ from the above expres- 
sions for the static masses, define the low-energy dynamical separation which characterizes 
the multicomponent quantum liquid p!5| , p!8|p3| . All this refers to the coherent part of the 



conductivity. The incoherent part is determined by the LWS II and non-LWS multiplets 

In this section we have presented the expression of the normal-ordered Hamiltonian and 
rapidity operators in the pseudoparticle basis of Tij. This confirms the consistency of the 
one-dimensional Landau-liquid theory which was shown to refer to this basis. The advantage 
of using the pseudoparticle basis is that the problem becomes perturbative, i.e. from the 
point of view of the pseudoparticle interactions it is possible to classify which scatterings 



are relevant. Consistent with earlier results [18] and our more detailed study in paper II, 



while at low energy and small momentum all pseudoparticle interactions are irrelevant, 
at low energy and large momentum only the two-pseudoparticle interactions are relevant. 
The same happens for low-energy excitations involving changes in the numbers N a . For 
simplicity and for further use in paper II, we have presented in this paper only the two first 
Hamiltonian terms of expression (39). However, Eqs. (30) — (32) contain full information 
about all the remaining terms of higher-scattering order. 
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In the companion paper II (see also [Tj|), the perturbative character of the present 
class of quantum liquids is used to study the symmetries and algebras at the critical point. 
These studies confirm that the Hamiltonian (39) is the correct starting point to construct a 
critical-point Hamiltonian. 



V. CONCLUDING REMARKS 



In the previous sections we have introduced and studied a new operator algebra fig 
describing the low-energy Hamiltonian eigenstates of integrable quantum liquids. Consid- 
ering explicitly the case of the Hubbard chain in a magnetic field and chemical potential, 
we have shown that the new operators create and annihilate the pseudoparticles of the one- 



dimensional Landau liquid p^Jlq] . Our algebraic approach permits an operator analysis of 
the BA solutions. Further, it allows us to construct a normal-ordered Hamiltonian which 
has in the pseudoparticle basis a universal form involving only k — 0, forward-scattering 
pseudoparticle interaction terms. Further, the perturbative character of the pseudoparticle 
basis implies that the corresponding two-pseudoparticle Landau / functions and forward- 
scattering amplitudes are finite, in contrast to the non-perturbative electronic representa- 
tion, in which the two-electron forward-scattering amplitudes and vertices diverge. To clarify 
further the context and implications of our results, we close with several remarks. 

First, concerning the relation to previous work on pseudoparticles [|l4| , [T5| . p!9i^0| , ^l| , we 
note that the pseudoparticle operator algebra motivates and justifies the validity of the 
Landau-liquid studies. These studies were based on explicit use of eigenenergies and eigen- 
functions, rather than the general algebraic approach. 

Second, although our explicit calculations were presented for the Hubbard model, which 
in the general terminology corresponds to a "two-component integrable quantum fluid" , our 



results in fact apply to the class of BA solvable models [p| pl3| , |T8f , which in the largest sector of 
the parameter space have symmetry [[/(l)] 1 ' ||18|| . As remarked previously, in this sector, the 
model has v independent branches of gapless elementary excitations. Therefore, when v > 1 
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these models are usually referred as multicomponent integrable quantum liquids |p| |U3| , |T8 . 
Note in this sector that v represents the number of independent conserved quantum numbers 
and also the number of external fields (magnetic field, chemical potential, etc.) which 
are "conjugate" to these conserved quantum numbers in the statistical mechanics sense. 
Although some of these systems are described by continuum models 0, most of them are 
defined in an one- dimensional chain with N a sites, j = 1, N a , refer to interacting fermionic 
particles, and, in particular, to interacting electrons, as in the case of the Hubbard model 
studied in this paper. 

More precisely, the class of solvable quantum liquids to which our results apply is defined 
by the fact that the v x v dressed- charge matrix [|3||13|,|l7^|18| is determined exclusively by the 



R matrix associated with the solution of the Yang-Baxter equation P,T8[- 111 this regard, we 



note that our study here, when supplemented by the results in II and Ref. shows that 
the present operator algebra is consistent with the finite-size-correction results of Frahm and 
Korepin who found that the complete critical theory of v multicomponent integrable 
systems is a direct product of v Virasoro algebras. 

Third, the pseudoparticle algebra allows us to see immediately the origin of the "univer- 
sal" character of these one- dimensional integrable quantum liquids, which can be understood 
as a straightforward generalization to pseudoparticles of Wilson's renormalization group ar- 



guments ||31|| : the pseudo- Fermi points of the pseudoparticles replace the particle Fermi 
points, and close to the pseudo- Fermi points only few types of two-pseudoparticle scattering 
processes are relevant for the low-energy physics. 

Fourth, concerning the very important problem of the relation of pseudoparticles to the 
original electrons, we should make several remarks. As we have stressed throughout, the 
essential simplifying feature of the pseudoparticle basis is that the quantum problem be- 
comes perturbative there, in the sense that it is possible to identify which pseudoparticle 
scatterings are relevant for the low-energy physics. For instance, in the new basis the two- 
pseudoparticle forward-scattering vertices and amplitudes do not diverge and one can define 
a many-pseudoparticle perturbation theory in which the non-interacting ground state is the 
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exact ground state of the many-electron problem. This is directly analogous to Landau's 



Fermi-liquid theory, in which the interacting ground state is used as reference state p6| , |37 
As is well known, the present class of quantum liquids are not Fermi liquids because the 
single-particle spectral function has no quasiparticle peak (i.e., Zp = 0), but they do have 
what was previously called a " Landau-liquid" character [Tj] , |n^|rj|gD| , gTJ : in the pseudoparti- 
cle basis the low-energy physics is fully controlled by two-pseudoparticle forward scattering, 
and writing the Hamiltonian and other operators in the pseudoparticle basis allows the use 
of perturbation-theory techniques, as we have shown in this paper and will exploit in II. 

The non-perturbative character of the usual electronic basis - e.g., the divergences of 
the two-electron forward-scattering vertices and amplitudes - is reflected in the complex 
and exotic properties of the electron-pseudoparticle operator transformation, which we have 
studiously avoided in this paper. The construction of this energy, momentum, and param- 
eter dependent nonlinear transformation is in general very involved. For instance, we know 
that at low energies and large momenta it maps one-pair electron operators onto multipair 



pseudoparticle operators, and vice versa [22" ] . Further, both the colors a and pseudomomen- 
tum q of the pseudoparticles are not, for general values of the parameters, simply related to 
the usual electronic quantum numbers (i.e., charge, spin, up spin, down spin, and momen- 



tum) |2!|,[34]]. We shall return to this problem elsewhere |34| and solve it for the generators 
which transform the exact ground state onto the low-energy and small momentum Hamil- 
tonian eigenstates, which we are able to write both in terms of pseudoparticle and electron 
operators. Importantly, as we have demonstrated in our calculations here, since the BA 
solution refers at low energy to the pseudoparticle operator representation, we can study 
the quantum problem in the corresponding basis without explicit knowledge of the complex 
electron-pseudoparticle transformation. 

Fifth, concerning applications of our ideas to other theoretical models and to problems in 
real materials, we note that although our present analysis refers to a class of one-dimensional 
integrable quantum liquids only, it has been argued that similar non-Fermi liquid behavior 
occurs in some sectors of parameter space of two-dimensional interacting quantum liquids 
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P J41| . |4^j43] . |44}1 . If this is so, techniques similar to ours may apply. Further, it is of consid- 
erable interest to determine the extent to which the algebraic structure and pseudoparticle 
perturbation theory remain valid in systems which are not integrable but which behave as 
Luttinger (or Landau) liquids (about the relation between Luttinger and Landau liquids 
see Sec. VIII of Ref. ||15| ); an example of such a system is the one-dimensional extended 
Hubbard model. In terms of real materials, it is known that the present one-dimensional 
quantum systems provide useful (albeit idealized) descriptions of the physics of quasi-one- 
dimensional solids p2| , ^5|,|46i . At low frequency the pseudoparticles are the transport carriers 
of the quantum liquid and couple to external fields |I^,[l5| . These couplings determine the 
exotic instabilities observed in quasi-one- dimensional synthetic metals ||22|| . 

Finally, we turn to directions for further research. In addition some specific calculations 
mentioned above, our immediate aims are to apply the concepts and techniques developed 
here to establish that the pseudoparticle perturbation theory leads to the correct effective 
Hamiltonian, which is obtained from the universal expression (39) — (40) of the quantum- 
liquid normal-ordered Hamiltonian in the pseudoparticle basis, and to study the Virasoro 



algebras [ |TB| , |i"8fl of conformal field theory. The initial results of these studies are described 
in II. The construction of a generalized Landau-liquid theory referring to the Hilbert space 
spanned by both the LWS I and LWS II is in progress. 
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APPENDIX A: NORMAL-ORDERED OPERATOR EXPRESSIONS 



Following the discussion of Sec. IV, the perturbative character of the system implies 
the equivalence between expanding in the pseudoparticle scattering order and/or in the 
pseudomomentum deviations (28). This, together with Eqs. (20) and (26), reveals that the 
problem of calculating the normal-ordered operator expansions (37) and (39) is equivalent 
to the Landau-energy functional studies of Refs. jl5|,^0l |2~lj . Since in these papers the full 
expression of the second-order function Q$(q) (which corresponds to the operator Q$(q)) 
and a number of other functions were not presented, we give in this Appendix a short 
description of the calculation of the normal-ordered operator expansions (37) and (39). We 
focus our brief study on the points which were omitted in Refs. [ I~5|,p0~ , pl| . 



We start by the evaluation of the first-order and second-order terms of the operator (37). 
For simplicity, we consider here pseudomomentum deviations and eigenvalues. Equation 
(26) then allows the straightforward calculation of the corresponding operator expressions. 

In the thermodynamic limit, Eqs. (31) and (32) lead to the following equations 

m - 1 /f; N M) t^(sw) - *2*!s>.) = „ (ad 



and 



1 r f q( Fc „ T , ls 1 f „, , sin K(q') 



'If, 

a (+) I 

q Fs _ , , N i/l 



, N s (q')tan-H-(S(q)-S(q')))] = q, (A2) 

iFs Z 

where we have replaced the operators by the corresponding eigenvalues and the summations 
by integrations and u = U/At. The eigenvalue form of Eq. (35) is 

5R a (q) = R° a (q + 5Q a (q)) - R° a (q) , (A3) 

where 5R a (q) and SQ a (q) are the eigenvalues of the operators : R a {q) : and : Q a {q) 
respectively. From Eq. (37) SQ a (q) can be written as 



5Q a (q) = Q { * ) (q) + Q ( a ) (q) + -, (A4) 
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where Q$(q) is the eigenvalue of the operator Q^(q). Expanding the 5R a (q) expression 
(A3) we find 

oo 

tfa(?)=E4 l) (?), (A5) 
i=0 

(and 5R a (q) = EZi R$(q)) where the 

zero-order ground-state functions (q) are defined 
by Eqs. (Al) and (A2) with the distributions N a (q) given by Eq. (24). From the resulting 
equations we can obtain all derivatives of the ground-state functions R^(q) with respect to 
q. The terms of the RHS of Eq. (A5) involve these derivatives. For instance, the first-order 
and second-order terms read 

R£\q) = d -^Q^(q), (A6) 

and 

*< 2) (9) = d -^Q i t>M + i^f^'to)] 2 , (A7) 

respectively, and involve the first and second derivatives. Let us use the notation of Eqs. 
(Al) and (A2), i.e., 

R c (q) = K(q) R s (q) = S(q), (A8) 

and 

Rf{q) = K^{q) Rf \q) = S®(q) . (A9) 

From Eqs. (Al) and (A2) (with N a (q') given by the ground-state distribution (24)) we 

find that the first derivatives dRa dq ^ can be expressed in terms of the functions R$(q) as 
follows 

dK^(q) 1 /A N 

, W = ?+j , A10 

dq 

and 



1 1 cosX(0)( g ) rq^J d , 1 

J q ( -J ^ l+[5(0)( g ')-^^2i] 
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^ = ^ ^ ■ (AH) 

dq 1 f ? Fc . , 1 L ( q Fs j n , 1 

- J q£ q 1+[S (0) {q) -^lm^ 2. J q ( F -) a( l l+[i(5(0)( 9 )-5W( g '))] 2 

From the same equations we find that the second derivatives d R J* q2 ^ can be expressed 
in terms of the functions R^(q) and its first derivatives (A10) and (All) and read 



<PK<®{q) _ { dK^ (q) ]3f l sin K^(q) fi 



= [— TT^n- - V 7 <-\ dq' 



(+) 



and 



dq 2 dq L vr u 7^ j + [£((%') - EEi^M ]2 

_ 2 cosing) [g (Q) (g0 _^M " 

vr [ « J ^ {1 + [5(0)^) - ™*<°>fa)]2 } 2*' lAi ' j 



rfg 2 rfg 7ri^- c ) {1 + - sin ^ 0) (g0 ]2} 2 

- -1 dq' ^\q)-S^q')\ ^ 

Let us introduce the functions Qc 1 ^) an d Q^H V ) sucn that 

Q«(?) = Q?(Kl%)) , Qf )((,) = Q?>(5<°)(9)) • (A14) 

Since we will often consider k and v integrations instead of q integrations, it is useful to 
define the following limiting values 

Q(±) = K {0) (q { pJ) , S W =5 (o) ((Z g) )> (A15) 

which refer to the pseudo-Fermi points (12) — (14). However, in order to be consistent with 
the order of the expressions below, these values are to be replaced by the corresponding 
leading-order terms ±Q and ±B, respectively, which read 

±Q = tf(°> (±q Fc ) , ±B = (±q Fs ) , (A16) 

where the pseudo- Fermi points ±qFa are defined in Eq. (15). Introducing both the distri- 
butions 

N a (q) = N° a (q) + 5N a (q), (A17) 
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and the first-order functions (A6) into Eqs. (Al) and (A2), we find after expanding to first 
order that the functions (q) and Qy' (q) are given by Eq. (A14) with Q^(k) and Q^(v) 
defined by the following system of coupled integral equations 

Q?\k) = 1 [f] dq>8N s {q>) t,n\S^\q') - — ) + - T dv'—^pQ^- , (A18) 

7T Jql > U 7T J-B 1 + \v' — SiSJIU 



and 



7T Jq)T> U 

+ - / dg'5iV s tan" 1 - u - 5 (q')}) + - / dfc' 



_|_ _ sin fc' j 2 

~ 2^ J-b l + [\{v-v')f ' (A19) 
The use of Eq. (A18) in Eq. (A19) allows the introduction of a single integral equation for 
Q^'(v). Combining the obtained equations with Eqs. (A14) and the phase-shift equations 



(23) — (38) of Ref. |2T| leads to expressions (38). 

In order to evaluate the second-order functions Q^\q) and Qy{q) of the RHS of Eq. 
(A4) we introduce the functions (A6) and (A7) and distributions (A17) in Eqs. (Al) and 
(A2). Expanding to second order we find after some algebra 

Q ( ^{q) = Q { ^\<l) + \j q m\<l)} 2 }, (A20) 

where 

Q^\q) = Q^\R^\q)), (A21) 

and the functions Q^*^(k) and Qf'*\v) are defined by the following system of coupled 
integral equations 

1 r* +) ^rt^dsww QW( q ') 



Qf*\k) = - dq'8N s (q 



nJqP ^ 1 "'^ s "' 1, dq' 1 + [SW{q') - 

1 r dS&>{q) IQ^(jqFs)} 2 
2n [ dq llq - qF ° + 

1 r» 

IT J-B 1 + \V' — 



+ r/-^ , , ,'„ _ l J»„ . (A22) 
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and 

,(+) 



^ 1 } TrAw aq0Iyc ^> dq > u 1 ] [t , sin^(°)( q >) ]2 

1 /^waw^m Q " )(g,) 



2tt ./,<-> 4 sv *' dg' l + fKv-W))] 

i r ^ (0) (g) n cosg [g«o^e)] a 

^ JI<?=9F C „, 2^ 



2tt l dq uq ~ qtc u ^ 1 + [ v - J -^} 2 
1 r dSW(q) [QP(jq F s)} 2 

1 [Q ,cos fc' Q^(A^) 



7T J-Q U l + \V — 



sinfc' 12 



2tt J- 



dv \^( w ' (A23) 

-B l + lziv-V 1 )] 2 

respectively. Note that the free terms involve the first-order functions. 

In order to derive the first-order and second-order Hamiltonian terms of Eqs. (40) and 
(41), we again consider eigenvalues and deviations. We introduce in the Hamiltonian expres- 
sion (30) Eqs. (A6), (A7), and (A17) and expand the obtained expression. Besides simpler 
terms, the first-order expression includes the term fSn dk[2t sin k]Q^ (k). After some algebra 
we find the first-order term of Eq. (40) with the pseudoparticle bands given by expressions 
(19) — (21) of Ref. PD| . Expanding (30) to second-order leads to the following energy 



E 2 = ^{ ff" dq'5N^)^^[2tsmK^(q')]Q^(q') 



+ 



dq' 

1 dK^\q) 



[tsinQ] E [Q { c\mfc)} 2 + I dk[2t sin k}Q^*\k)} , (A24) 

.•ZTi J-Q 



2n l dq Jl <™ 

such that \r) z , S z ) = E 2 \ri z , S z ). Inserting the suitable functions in the RHS of Eq. (A24), 
performing some integrations by using symmetry properties of the integral equations (A22) 
and (A23), and replacing deviations by pseudomomentum normal-ordered operators (25) we 
find after some algebra expression (41). This can be rewritten in terms of the / functions 
(42) as given in the RHS of Eq. (40). 

* Permanent address, Department of Physics, University of Evora, Apartado 94, P-7001 
Evora codex, Portugal. 
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